Ortho-to-Para Ratio Studies of Shocked H2 Gas in the Two 
Supernova Remnants IC 443 and HB 21 

Jong-Ho Shinn 1 , Ho-Gyu Lee 2 , Dae-Sik Moon 3 
jhshinnOkasi .re .kr , hgleeOastron. s .u-tokyo .ac.jp, moon@astro .utoronto . ca 

ABSTRACT 

We present near-infrared (2.5 — 5.0 //m) spectral studies of shocked H 2 gas in 
the two supernova remnants IC 443 and HB 21, which are well known for their 
interactions with nearby molecular clouds. The observations were performed with 
Infrared Camera (IRC) aboard the AKARI satellite. At the energy range 7000 
K < E(v, J) < 20000 K, the shocked H2 gas in IC 443 shows an ortho-to-para 
ratio (OPR) of 2.4±g;|, which is significantly lower than the equilibrium value 3, 
suggesting the existence of non-equilibrium OPR. The shocked gas in HB 21 also 
indicates a potential non-equilibrium OPR in the range of 1.8 — 2.0. The level 
populations are well described by the power-law thermal admixture model with 
a single OPR, where the temperature integration range is 1000 — 4000 K. We 
conclude that the obtained non-equilibrium OPR probably originates from the 
reformed H2 gas of dissociative J-shocks, considering several factors such as the 
shock combination requirement, the line ratios, and the possibility that H2 gas 
can form on grains with a non-equilibrium OPR. We also investigate C-shocks 
and partially-dissociative J-shocks for the origin of the non-equilibrium OPR. 
However, we find that they are incompatible with the observed ionic emission 
lines for which dissociative J-shocks are required to explain. The difference in the 
collision energy of H atoms on grain surfaces would make the observed difference 
between the OPRs of IC 443 and HB 21, if dissociative J-shocks are responsible 
for the H 2 emission. Our study suggests that dissociative J-shocks can make 
shocked H 2 gas with a non-equilibrium OPR. 
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1. Introduction 



1989; Kaufman fc Neufeld 



Infrared H2 line emission is a useful tool for studying shock-cloud interactions. H2 is one 
of basic coolants at the post-sho ck region (cf. Neufeld Sz Dalgarnolll989 ; Hollenbach fc McKee 



1996), and its infrared line emission usually suffers relatively small 



extinction effects (cf. lDraindl2003l ). The shock-excited H 2 line emission has been widely ob- 



served from diverse t argets, such as pro t ostellar Jets a nd outflows, supernova remnants, and 



interacting galaxies (IHabart et al.ll2005t IOmontll2007t and references therein). Ground tele 



scopes usually observe high excitation emission lines (E(v, J) > 7000 K) in the near-infrared 
wave bands, while space satellites observe low excitation emission lines {E(v, J) < 7000 K) 
in mid-infrared. The high excitation lines show the excitation temperature (T ex ) of a few 
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2001: Giannini et al. 2006: Caratti o Garatti et al. 2008 



H 2 line emission manifests one important diagnostic parameter of H 2 gas: the ortho-to- 
para ratio (OPR). First we give a brief description of ortho- and para-H 2 . An H 2 molecule 
consists of two identical proton nuclei, i.e., two identical fermions; hence, for a given nu- 
clear spin combination, a specific rotational quantum n umber (space wave-function) is only 
perm itted according to the Pauli exclusion principle (cf. Field et al.lll966t IShull fc Beckwith 
19821 ). If the nuclear spins are summed to zero (/ = 0, singlet), it is called as para-H 2 and 
only even rotational quantum number (J) is allowed. On the other hand, if the nuclear spins 
are summed to one (1=1, triplet), it is called as ortho-H 2 and only odd- J is allowed. The 
statistical weight (gj) for a level (v, J) is a multiplication of the nuclear spin degeneracy 
and the rotational degeneracy, gj = (21 + l)(2J + 1); i.e., gj = 3(2J + 1) for ortho-H 2 and 
gj = (2J + 1) for para-H 2 . 

The OPR is the ratio between the abundance of ortho- and para-H 2 in all v, J levels. 
In general, lower E(v, J) levels are much mor e populated than higher E(v, J) levels (cf. 
Timmermann et al.lll996t iRosenthal et al.ll2000l ). Therefore, the OPR is mainly determined 
by the OPR at low E(v, J). This means that the OPR can be different from the abundance 
ratio of ortho- and para-H 2 in high E(v, J). In this paper, we make a note on the energy 
level E(v, J) when we mention the OPR at a specific E(v, J) range; otherwise, the OPR 
means the abundance ratio of ortho- and para-H 2 summed over all v, J. We also note that 
the OPR of H 2 in local thermodyn amic equilibrium (LTE) is equal to 3 at T > 300 K, 
decreasing from 3 at T < 300 K (cf. burton et aDll992l : Isternberg fc Neufeld! ^9991 ) . When 
the OPR is different from this LTE OPR, we call it as a non-equilibrium OPR. 



The OPR of H 2 is useful in estimating thermal history of the gas ( jNeufeld et al.ll2006l ). 
Under typical conditions of molecular clouds, the collisional interconversion between ortho- 
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and para-H2 is much slower than the collisional transitions between the same-type H2 (ortho- 
to-ortho and para-to-para). Therefore, the OPR have the information on how much inter- 
conversion collisions the H 2 gas experiences and what the initial OPR of the H 2 gas is. 
For shocked H 2 g as, th e possibility of the non-equilibrium OPR was briefl y discussed by 
Chang fc Martin dl99lh. and later this question was seriously addressed by iTimmermannl 



(119981 ) and IWilgenbus et al.l ( 120001 ) . The last two studies predicted the existence of H 2 gas 
with the non- equilibrium OPR in C- shocks or partially-dissociative J-shocks (for the types 
of shocks, see lDraine fc McKedll993l ). 



Hitherto observations have revealed shocked H 2 gas of non-equilibrium O PRs from mid- 
infra r ed observations, which trace low energy levels of E(y , J) 5, 7000 K (e.g., iNeufeld et al. 
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Giannini et al.l l2006t ICaratti o Garatti et al.l 120061 ) . Here we report two cases suggesting 
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the shocked H 2 gas with non-equilibrium OPRs at such high levels, which is revealed from 
near-infrared spectral observations with the AKARI satellite. The targets are two super- 
nova remnants IC 443 (G189. 1+3.0) and HB 21 (G89.0+4.7), which are w ell-known for their 



interaction s with nearby molecular clouds (for more about the targets, see lShinn et al.ll2009 



20101 . l201ll ). We present new HB 21 data here, while we used the previously pr esented data 



for IC 443 (clump B; i.e., IC 443B), which had been left for future analyses (jShinn et al. 



20111 ). The results are discussed with respect to different shock types, and we conclude that 
the obtained non-equilibrium OPR may originate from the reformed H 2 gas of dissociative 
J-shocks. 



Observations and Data Reduction 



The AKARI near-infrared spectroscopic observations were performed on two regions 
of HB 21, where the shock-cloud interactions are under way. The two regions, "HB 21N" 
and "HB 21S", ar e listed in Table ED The " HB 2 IN" and "H B 21S" correspond to the 
regions "N2front" (jShinn et al.ll2009l ) and "SI" ( jShinn et al.ll2010l ). respectively. Two nearby 
background positions ("HB 21N BG" and "HB 21S BG," Tabled]) were also observed in order 
to inspect if any background H 2 line emission exists. 



The observations were carried out with the Infrared Camera (IRC. lQnaka et al.ll2007l ) of 
the AKARI satellite (IMurakami et al.ll2007l ). on 2008 June 1-3 during the post-helium phase. 
All the observations were performed with the observation mode IRCZ4, which is designed 
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for general spectroscopic observations. The IRCZ4 mode acquires the spectra from multiple 
exposures without dithering, taking one reference image in the middle of the exposures 
( lOnaka et al.l 120091 ). Based on this refere nce image, we correcte d the observed astrometric 
positions employing the 2MASS catalog (ISkrutskie et al.l 120061 ); the resulting astrometric 
accuracy is smaller than 1.4". We used the 5" x 48" slit and the grism, whose spectral 
resolution and wavelength coverage are AA ~ 0.03 /im and 2.5 — 5.0 jum, respectively. Figure 
[U shows the slit positions over the H 2 v = 1 — > S'(l) 2.12 /im image of HB 21. The H 2 
2.12 /im image was acquired with the Wide-field InfraRed Camera (WIRCAM) of Canada- 
France-Hawaii Telescope (CFHT) on 2007 August, and the data reduction was carried out 
employing the SIMPLE package (jWang et al.l 120101 ). 

We fi rst applied the re duction pipeline for the AKARI post-Helium data (version 
20110301. lOnaka et al.ll2009l ). The column pull-down effect was corrected and the hot pixels 
were masked out. No smoothing and tilt-correction were applied to the two dimensional 
spectral image, while running the pipeline. Then, we extracted a spectrum for each obser- 
vation ID (cf. Table [T|) by averaging the pixel values along the slit direction. Again these 
multiple spectra for each pointing position were averaged together, and the results are pre- 
sented in Figure [2J Some spectra show an alternating noise pattern along the wavelength 
bin, but they do not affect the measurement of line intensity we intend. The calibration 
systematic error of 10 % (private communication with the AKARI helpdesk) was squarely 
summed to the line intensity error (Table |2]). 



3. Analysis and Results 

3.1. Line Intensity and Level Population 

Figure |2] shows the AKARI spectra of HB 21. A few evident emission lines are seen 
from the spectra of "HB 21N" and "HB 21S." On the other hand, no emission line is seen 
from the spectra of the background positions, "HB 21N BG" and "HB 21S BG"; hence, the 
background contribution of line emission is negligible. We identified all th e detected emissio n 



lines are of H 2 , comparing the line positions to those detected in IC 443 ( IShinn et al.l 120111 ). 



The line intensities were measured by fitting a Gaussian profile plus a synthetic contin- 
uum to the emission line spectrum. The synthetic continuum was made by median-smoothing 
the observed spectrum with a box kernel. The kernel width is wide enough to smooth the 
line features out; it is about 0.19 /xm. This continuum is displayed as a dashed line in Fig- 
ures [3] and |H We set the full-width-at-half-maximum of the Gaussian profile as the spectral 
resolution of IRC, AA = 0.03 /im. When two lines are overlapped, we fitted both lines simul- 
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taneously. The fitting results are shown in Figures |3] and HI and the measured line intensities 
are listed in Table [2j 

In order to derive the level population of the shocked H 2 gas, we first performed the 
reddening correction to the measured line intensities . We employe d the extinction curve of 
"Milky Way, R v = 3.1" (TWeingartner & Drainel I200I lDrainell2003[ ). and adopted N U =N(R 
I)+2N( ~H.<>)= 3.5 x 10 24 cm -2 derived from the X-ray absorption towards the central region of 
HB 21 (ILee et al.ll200ll ). Then, these reddening-corrected line intensities were used to derive 
the level population, assuming that the H 2 emission is optically thin. This assumption is 
appropriate since the H 2 emiss ion lines become optically thick when iV(H 2 ) > 10 24 cm -2 , 
as noted in IShinn et al.l (1201 ll ) . If th ere was such a hig h iV(H 2 ) > 10 24 cm -2 toward HB 
21, the AKARI mid-infrared images (IShinn et al.l 120091 ) should have shown absorbed fea- 
tures against backgroun d radiation, because the iV(H 2 ) exceeds that of infrared dark clouds 
faerem fe Tafallall2007h . 



We used the following equation to derive the extinction-corrected level population of 
the shocked H 2 gas. 

at < n 4?rA I rc (v, J -+ v', J') 

Nrc ^ J) " A(v, J -+ V, JO ' (1) 

I rc (v, J — > v', J') and A(v, J — > v', J') are the reddening corrected line intensity and the 
Einstein-A radiative transition probability, respectively. A, h, and c are the transition wave- 
length, Planck constant, and light speed, respectively. We employed t he molecular physica l 
constants tabulated in the simulation code CLOUDY (version C08.00; iFerland et al.lll998l ). 
The derived column densities are listed in Table |3j Their population diagrams are shown in 
Figure note that the LTE population appears as a straight line in this diagram. 

The populations of (v, J) = (0,11), (0,12), (0,13) show a zigzag pattern. The OPRs 
estimated with the single T ex assumption are 1.8±0.4 and 1.7±0.3 for HB 21N and HB 21S, 
respectively (T ex ~ 2000 K). These OPRs ar e non-equilibrium OPRs (cf. section [1]). Note 
that IC 443B also showed the zigzag pattern (IShinn et al.l 120 111 ); the OPRs are estimated to 
be 2.4 ± 0.4 and 2.1 ± 0.3 for (v, J) = (0, 11-13) and (1,1-3), respectively (T ex ~ 1000 - 2000 
K). We here underline that these (v, J) levels all reside at high energy, E(v, J)> 7000 K. 

In order to ascertain that the zigzag pattern seen from the levels (v, J) = (0, 11-13) reflects 
the intrinsic OPR < 3.0, we examined as follows. We checked the flux calibration procedure 
including the data reduction pipeline, and found no indication of systematic features that 
can mimic the zigzag pattern. We also checked if any blended line exists for v = — S(10) 
4.41 fim, which corresponds to the level (v , J) = (0,12). Ther e is no known strong atomic or 
ionic lines from shocks around 4.41 fim (cf. lAUen et al.ll2008l ). The line blending between H 2 
emission lines is dealt with in section 13.2.11 
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3.2. Power-law Thermal Admixture Model of Shocked H 2 Gas with a 

Non-equilibrium Ortho-to-Para Ratio 

The power-law the rmal admixture m odel has been used to describe the level popula- 
tion of shocke d gas rtShinn et al.ll201ll . and references therein). We employed the same 

( 201lh . and slightly modified it for the reproduction of the non- 



model used in Shinn et al. 



equilibrium OPR, i.e., the zigzag pattern in the population diagram. The relative pop- 
ulation between ortho- and para-levels is controlled by scaling the rate coefficients for the 
reactive collisions with H, which mediate the ortho-to-para and para-to-ortho transitions (cf. 
Le Bourlot et al.lll999l ). The final level population is obtained by integrating the following 
equation, 



dN 



aT~ b dT 

N(# 2 ;T > 1000K)(6- 1) 

rpl— b _ rpx—f) 
min max 



(2) 
(3) 



where a and b are constants and T min = 1000 K and T max = 4000 K. The statistical 
equilibrium was assumed and the collisional partners are H2, He, and H. We assumed 
n(He)= 0.2n(H 2 ), and the abundance ratio of H to H 2 was set as a free parameter Xh- 



Based on the previous observations of IShinn et al.l (|20111). we adopted T , 



1000 K 



i nstea d of 100 K unlike previous applications (e.g., IShinn et al.l 120091 . 120101 ) . IShinn et al. 
( 1201 if ) found that two different densities are required to describe the level population of 
shocked H 2 gas in the range of < E(v, J) < 25, 000 K: one with n(H 2 ) of ~ io 2 - 8-3 - 8 cm -3 
for E(v, J) < 7000 K and the other of ~ io 5 - 4 " 5 - 8 cm" 3 for E(v, J) > 7000 K. From these 
results, they recognized that only T > 1000 K gas is crucial to describe the level population 
of E(v, J) > 7000 K. Therefore, we adopted T min = 1000 K for the model calculation, since 
our data points spread over E(v, J) > 7000 K (Figure E]). 

We then applied the model to the reddening-corrected level population of shocked H 2 
gas in IC 443B, HB 21N, and HB 21S. The fitting results are shown in Figure [61 and the 
corresponding model parameters are listed in Table HI In the case of IC 443B, we first 
tried to fit the data without any fixed parameters, but failed to have inversely diverging 
n(H 2 ) and Xh', ^(H 2 ) increases over 10 7 cm -3 and Xh decreases below —3.0. These two 
parameters seem to be not well-constrained by the given data and error. Since such a high 
ri(H 2 ) is unreasonable, we fixed the parameter Xh- Two different fittings were performed 
with Xh = —1.7 an d Xh = —2. 5. Xh = —1.7 is chosen since it was obtained from the 



analysis of IC 443C ( IShinn et al.l 1201 ll ). Xh = —2.5 is chosen since IC 443B seems to be 
experiencing a weaker shock than IC 443C, where Xh = —1.7 was obtaine d. The intensit y 
ratio /(CO; J = 3 - 2)//(CO; J = 2 - 1) is lower at IC 443B than IC 443C JXu et al.lboilh . 
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and the "W-shaped" H 2 morphology of IC 443 suggests that IC 443B may be denser than 
IC 443C (see Fig. 1 in lshinn et allboilh . In the case of HB 21N and HB 21S, we fitted the 
data fixing both n(H 2 ) and X H , because the data points are too few. We just show the level 
population can be described by the power-law thermal admixture model. We adopted n(H 2 ) 
= 10 5 5 cm " 3 and Xh = —2.0 , which are thought to be typical values for shocked H2 gas of 



IC 443 (cf. IShinn et al.l 120111 ) . Figure [7] shows the \ 2 contour of the model parameters for 
IC 443B. We did not perform the x 2 contour scan for HB 21N and HB 21S, since the data 
points are too few, compared to the number of model parameters. 



IC 443B s hows n(H a ) = 10 



5.8-6.5 



10 



cm 3 somewhat higher than that of IC 443C, n.(H 2 ) 



4 cm 3 (IShinn et al.l 1201 ll ). This density difference may be related with the "W- 



shaped" H 2 morphology mentioned above. iV(H 2 ) of IC 443B is 10 cm ; here, 

A r (H 2 ) means iV(H 2 ; T > 1000 K), i.e., the column density of H 2 gas whose temperature 
is higher than 1000 K. IC 443B shows b = 3.9 - 4.5 larger than that of IC 443C, b = 1.6. 
This difference is natural because b of IC 443C is mainly determined by the data of E(v, J) 
< 7000 K, while it is not for b of IC 443B. The obtained OPRs are 2.1 - 2.2. Interestingly, 
the OPR of IC 443B shows a smaller variation than the other parameters under the different 
adoptions of (cf. Table H]). These OPRs are similar to those estimated from the levels 
(v, J) = (0, 11-13) and (1,1-3) assuming a single T ex : 2.4 ± 0.4 and 2.1 ± 0.3. The \ 2 contour 
(Figure[7]) shows that there is a weak anti-correlation between n(H 2 ) and b, but no correlation 
between iV(H 2 ) and the OPR. Figure [7J definitely shows the non-equilibrium OPR even in 
95% confidence level. 



HB 21N and HB 21S show iV(H 2 ) of 10 



, 4 cm 2 



and 10 cm , respectively. This is 



about an order of magnitude smaller than that of IC 443C, which explains the consistent 
difference in the line intensities (Figure |2j IShinn et al.l 120111 ) . HB 21N and HB 21S show 
b values of 2.6 and 3.3, smaller than IC 443B; this difference is hardly meaningful because 
the number of HB 21 data points is just half of IC 443B. The OPRs are 1.8 and 1.6 for HB 
21N and HB 21S, respectively. These values are similar to those estimated from the levels 
(v, J)=(0,ll-13) assuming a single T ex : 1.8 ± 0.4 and 1.7 ± 0.3. 



3.2.1. Spectral Fitting to the Observed Spectrum 



In order to investigate the OPR of shocked H 2 gas in IC 443B using all possible obser- 
vational constraints, we additionally performed spectral fitting to the observed AKARI-IRC 
spectrum of IC 443B. The model spectrum is made by convolving the modeled intensity 
with a Gaussian kernel whose full-width-at-half- maximum is 0.03 fim. The extinction effect 
is included using Ay = 13.5 (IShinn et al.l 1201 lh and the extinction curve of "Milky Way, 
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Rv = 3.1" (jWeingartner fe Draindl200ll ; lDraindl2003l ). We used a synthetic continuum that 
is made by median-smoothing the observed spectrum with a 0.29 /iin box kernel, and the 
continuum is fixed during the fitting. The calibration systematic error of 10 % (cf. section 
[2]) is quadratically included in the error of spectrum. The wavelength shift is included as 
a free parameter, and found to be ~ +0.008 /im for the best fit. The fitting results are 
shown in Figure [HJ We fitted the spectrum fixing Xh = —1.7 and Xh = —2.5 as in the level 
population fitting (Fig. [6]). The best fit model parameters and the y 2 contours are shown in 
Table \5\ and Figure M 

The OPR is found to be 2.4$% which is a little higher than the ones from the level 
population fitting (Table HJ). However, the OPR is still less than 3.0 with a 95% significance 
limit (Fig. [H]). The index b also increases a little, while iV(H 2 ) and n(H 2 ) are almost the 
same. The difference between two kinds of fitting is likely induced by the constraints from 
the weak and blended H2 emission lines, which are ignored in the level population fitting 
(Fig. [6]). Such contributions can be seen in Table [6j for example, the three lines around 4.41 
/an, v = - S(10), v = 1 - 1 S(ll), and v = 2 - 1 0(8). 

For the HB 21 data, we could not get the range of OPR, because there are only a few 
emission lines strong enough to constraint the model parameters. We estimate the OPRs of 
HB 21N and HB 21S with a single T ex assumption using the 4.2, 4.4, 4.7 /im emission lines 
in section [341 They are found to be 1.7 — 1.8. If the contribution of v — 1 — 1 S(ll) 4.42 
fim to the 4.4 /im line feature and the contribution of v = 1 — 0(8) 4.16 /tm to the 4.2 
am. line feature are respectively around 13% and 10% as in IC 443B (cf. Table E]), then the 
OPRs increase to 1.8 — 2.0. However, this estimation is still tentative. More near-infrared 
H 2 emission line data are required to assess the OPR of shocked H 2 gas in HB 21. 



4. Discussion 

4.1. Origin of the Non-equilibrium Ortho-to-Para Ratio Obtained from E(v,J) 

> 7000 K 



We obtained the results suggesting shocked H 2 gas of non-equilibrium OPRs in IC 443B 
(OPR= 2.4t°j), HB 21N and HB 21S (OPR= 1.8-2.0, cf. section [3X1]). First, we examine 
the possibility th at this non-equilibr ium OPR originates from dissociative J-shocks, based 
on the results of Ishinn et all (boilh . They studied shocked H 2 gas in IC 443 (IC 443C 
and IC 443G), and found that the levels we here study (7000 < E(v, J) < 20000 K) are 
likely excited by dissociative J-shocks. Dissocia tive J-shocks mainly excite those levels at 
the H 2 reformation zone (cf . iFlower et al.l 120031 ) through two main excitation mechanisms: 
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the H2 formation pumping and the collisions with ambient gas ( IHollenbach fc McKedll989l ). 
Therefore, the non-equilibrium OPR most likely originate from the H 2 reformation zone, if 
dissociative J-shocks are responsible for the level excitation of 7000 K < E(v, J) < 20000 K. 

When the levels of 7000 K < E(v, J) < 200 00 K is dominantly populated, the tempera- 
ture is a few hundred kelvin (IFlower et al.ll2003l ). This temperature is too low to render the 
ortho-para interconversion. For ex ample, the number of para-to-ortho conversion is about 



0.6 for the typical J-shock seen in IFlower et al.l (120031). W e use the rate coefficient k po of 



8 x 10~ n exp(-3900/T) cm 3 s" 1 for H I (ISchofieldlll967f ). n(H I)= 10 6 cm' 3 , t = 100 yr, 
and T = 300 K. Therefore, the OPR upon reformation behind dissociative J-shocks would 
remain unchanged until the gas cools down as low as the preshock gas. This means that the 
non-equilibrium OPR should be acquired when H 2 is reformed, if dissociative J-shocks are 
responsible. 



H-2 reforms on grain surface s or through the gas ph ase reaction (Ivan de Hulstl 11948 
McDowelllll96ll ; iField et al.lll966l IShull fc Beckwithlll982l and references therein). Within 
the H 2 reformation zone of dissociative J-shocks, the grain surface reformation is domi- 
nant, because the gas phase reformation is ine fficient when the ionization fraction of shocked 
gas is below 2% (INeufeld &: Dalgarnd Il989l ). The OPR of reforme d H 2 was usually as- 
sumed to be the equilib rium OPR in dissociative J-shock models (IHollenbach fc McKee 



1989 



1987 



Flower et al.ll2003|). following the studies on the level population of newly formed 



H9 (IBlack fe Dalgarnolll976l ; iBieniek fc Dalgarnolll979t iBieniek 



1980; Black fc van Dishoeck 



Flower et al.l 120031 ). Thus, we need to more scrutinize previous theoretical and exper- 



imental studies relevant to the OPR of newly formed H 2 . 

As far as we are aware, there are two studies that pe rformed a quantum calculation for 



the H ? formation upon grains and mentioned the OPR (IMeijer et al.l l200ll ; iMorisset et al. 



20051 ). Among them, the results of IMorisset et al.l (120051 ) do not ex clude the pos s ibility that 
the newly-formed Ho has an OPR = 1.8 - 2.4 (cf. section ETU- iMeiier et all toOlh and 
Morisset et al.l (120051 ) respectively studied the Eley-Rideal (ER) and Langmuir-Hinshelwood 
(LH) formation mech anism^, among which the LH mechanism is more relevant to the inter- 
stellar environments (IPirronello et al.ll2000f ). The level population upon formation is biased 
within v > 5 for the LH mechanism, while it is scattered over all the considered vibrational 
levels (v = — 6) for the ER mechanism. They both calculated the OPR, and it is fluctuating 
within ~ 1.5 — 6.0 over the collision energy of 4 — 50 meV in the LH mechanism, while it 
is close to 3 over the initial translational energy of — 200 meV in the ER mechanism. We 
think that the OPR= 1.8 — 2.4 can be acquired through the LH formation mechanism, if the 



for the ER and LH mechanisms. 
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collision energy of H atoms on grain surfaces has an appropriate distribution. Note that the 
biased level population (v > 5) of the LH mechanism is redistributed over all levels through 
collisions with ambient postshock gas (i.e. ortho-to-ortho and para-to-para transitions). 

We found two experimental works tha t studied the OPR of newly-formed H2 on grains 
(jPerry fc Pried 120031 ; ICreighan et al.l 120061 ) . However, their results cannot tell if the newly- 
formed H 2 can have an OPR = 1.8 — 2.4 or not, because of the experimental contaminations. 
They measured the level population of newly-formed H 2 on a graphite surface, called "highly 
oriented pyrolytic graphite," at a base pressure of ~ 10~ 10 Torr (~ 1 9 cm" 3 K) much highe r 
than the typical pressure of interstellar clouds, 10 4 — 10 6 cm -3 K (jBergin fc Tafallal 120071 ). 



The obtained population is close to the LTE. However, it was known that the measured H 2 
populat ion had been cont aminated by the undissociated H 2 that originates from the H atom 
source (jlslam et al.ll2010l ). 



From the above discussion, we show that the non-equilibrium OPRs of 1.8 — 2.4 can 
originate from dissociative J-shocks. We here note that the dissociative J-shock interpre- 
tation is not incompatible with the observed line ratios between H 2 and H I. We compare 
the intensities of Bra 4.05 fim to those of v = — S(ll) 4.18 fim. Table UJ shows the 
expected and observed intensity of Bra. The expected Bra intensity is calculated from the 
observed v = — S(Tl) intensity (IShinn et al.l 1201 ll . Table [2]), employing the dissociative 
J-shock model of iHollenbach fc McKed (119891 ) . We ignore the effect of extinction, because 
the wavelengths of Bra and v — — S(ll) are similar. The adopted shock model parame- 
ters are n = 10 3 cm -3 , v s = 40 — 80 km s^ 1 . Note that the di ssociative J-shocks are likely 
propagating into interclump media of < 10 3 cm -3 , according to lShinn et al.l (120111 ). As seen 
from Table UJ, the observed intensities better agree with the model expectations, when the 
shock velocity is slower. The intensities ofi> = — 0S(11) are also compatible with the 
model values, considering the possibilities of the multiple shock surfaces, the shock surface 
viewed edge-on, the not-fully-covered shock surface because of unresolved substructures, etc. 

Now we examine other possibilities than the dissociative J-shocks. The shocked H 2 
gas of non-equilibrium OPR can be made by C-shocks or partially-dissociative J-shocks 
(jTimmermannl Il998l ; IWilgenbus et al.ll2000l ). However, thes e shocks are incap able of ema- 
nating the ionic lines observed in the shocked gas. When IShinn et al.l (120111 ) interpreted 
the level population of shocked H 2 gas in IC 443C and IC 443G (0 < E{v, J) < 25000 K), 
they concluded that a combination of C-shocks and dissociative J-shocks is needed. Under 
this interpretation, the dissociative J-shocks are mandatory to explain the observed ionic 
emission lines, like [Ne III 12.8 /zm and [C II] 158 /im. Such lines were also observed at 
IC 443B (jNoriega-Crespo et al. 2009 ). In HB 21, no detection of such ionic lines has been 
reported. However, the level population shape of shocked H 2 gas over K < E(v, J) < 20000 
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K is similar to that of IC 443C and IC 443G (Fig. EJ IShinn et al.ll2009l . bold l201lh — the 
power-law thermal admixture model with two density components is required to reproduce 
the level population. As in IC 443, dissociative J-shocks may be needed to explain the level 
population of shocked H 2 gas in HB 21. We thus think that the non-equilibrium OPR is 
unlikely to originate from either C-shocks or partially-dissociative J-shocks. 



4.2. Ortho-to-Para Ratios of Shocked H 2 Gas in IC 443 and HB 21 



The OPR obtained from the levels of 7000 K < E(v, J) < 20000 K is probably in non- 
equilibrium for IC 443B (section 13.2.11) . while it is in equilibrium for IC 443C and IC 443G 
(IShinn et al.ll201ll ). If we accept the conclusion from section |4TT| these OPRs are determined 
by the collision energy distribution of H atoms on grain surfaces. Therefore, the resultant 
OPRs means that the collision energy distribution is similar each other at IC 443C and IC 
443G, while it is different for IC 443B. It is uncertain how the collision energy distribution 
is determined. The property of grain surface and the environment where grains and H I gas 
are embedded may be included in the determining factors. We here note that IC 443B is 



located at the vertex of t 



(cf. Fig. 1 of IShinn et al 



than the others (cf. IShinn et al.ll201ll ) 



he "W -shaped" H 2 feature, while IC 443C and IC 443G are not 
20 1 1 ) . Thi s implies that the density at IC 443B may be higher 



The OPRs of HB 2 IN and HB 21 S are similar each other (Table H]). This indicates 
the collision energy distribution is similar at both regions, if we accept the conclusion from 
section 14.11 The OPRs of HB 2 1 are lower than those of IC 443. HB 21 shows a lower 
H 2 column density than IC 443 JNeufeld et alil2007l : Ishinn et al.lhood boiol . boilh . This 
implies a lower density for HB 21, if a similar extension of shocked H 2 gas along the line of 
sight is assumed for IC 443 and HB 21. 

If our inference for the preshock density of IC 443 and HB 21 is the case, the preshock 
density and the obtained OPR show neither a correlation nor an ant i- correlation. They 
are anti-correlated for IC 443B, IC 443C, and IC 443G, while they are correlated for IC 
443 and HB 21. This indicates that the obtained OPRs are not simply determined by the 
preshock density. More information on physical parameters is required to understand how 
the obtained OPRs are determined at high energy levels of 7000 K < E(v, J) < 20000 K. 
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5. Conclusions 

We present the near-infrared (2.5 — 5.0 /im) spectral observations of shocked H2 gas in 
the two SNRs IC 443 and HB 21. The observed regions are IC 443B, HB 21N, and HB 21S. 
At the energy range 7000 K < E(v, J) < 20000 K, IC 443B shows an OPR of 2.4i°J (< 3 in 
95% confidence limit), which suggests the existence of non-equilibrium OPR. HB 21N and HB 
21S also show indications of non-equilibrium OPRs in the range of 1.8 — 2.0. These may be 
the first objects observed with non-equilibrium OPR at such a high E(v, J) range. The level 
populations are well described by the power-law thermal admixture model (1000 <T< 4000 
K) with a single OPR, where the infinitesimal column density is proportional to the inverse 
power of the temperature. 



We conclude that the obtained non-equilibrium OPR probably ori ginates from t he H 



2 



reformation through dissociative J-shocks, based on several arguments. IShinn et al.l (120111 ) 
found that the H 2 levels of 7000 < E{v, J) < 20000 K are likely excited by dissociative 
J-shocks, where H 2 is dissociated and then reformed. Previous quantum calculations do not 
exclude the possibility that the reformed H 2 has an OPR of 1.8 — 2.4, when formed with an 
appropriate collision energy of H atoms on grain surfaces. The observed line ratio between 
H 2 -u = — 0S(11) and Bra is also compatible with the dissociative J-shock interpretation. 

According to our conclusion, the difference in the obtained OPRs is caused by the 
different collision energy of H atoms on grain surfaces. It is uncertain how the collision 
energy is determined, but the property of grain surface and the environment where grains 
and H I gas are embedded may be included in the determining factors. We check if any 
simple relation is held between the obtained OPR and preshock density for IC 443 and HB 
21, but it turns out unlikely. 

More observational studies are required to confirm the non-equilibrium OPR; to validate 
our explanations for the obtained non-equilibrium OPR; and to grasp how a non-equilibrium 
OPR of shocked H 2 gas is determined at the energy levels of 7000 K < E(v, J) < 20000 K. Our 
study suggests that dissociative J-shocks can make shocked H 2 gas with a non-equilibrium 
OPR. 



This work is based on observations with AKARI, a JAXA project with the participation 
of ESA. The authors thank all the members of the AKARI project. The authors also 
appreciate the referee's comments, which improve the original manuscript. J. H. S. is grateful 
to Bon-Chul Koo and Michael Burton for the discussion that clarified the ortho-to-para 
ratio and related modeling, and grateful to Ji Yeon Seok for the useful discussion on the 
data reduction. This research has made use of SAOImage DS9, developed by Smithsonian 
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Astrophysical Observatory (jJoye fe Mandelll2003l ) 
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RA (J2000) RA (J2000) 



Fig. 1. — The observed sli t positions for HB 2 IN (left) and HB 21 S (righ t). The two regions 
correspond to "N2front" Jshinn et alJl2009h and "SI" Jshinn et al.lboioh . respectively. The 
background images are H2 v = 1 — >• S(l) 2.12 /im images obtained with WIRCAM of 
CFHT telescope. The direction to the remnant's center is southward and northward in the 
left and right figures, respectively. 
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Fig. 2.— The AKARI IRC near-infrared spectra of HB 21N (top), HB 21S (bottom), and 
their backgrounds (cf. Fig. [1]). The vertical bars indicate the location of H 2 emission lines 
(cf. Table [2]). "?" indicates the location of doubtful H 2 emission line. For clarity, the error 
bars are omitted; the errors estimated from the pipeline are shown in Figures [3] and HI 
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Fig. 3. — Fitting for the H 2 emission lines observed toward HB 21N. The solid and dashed 
lines indicate the continuum+line and continuum, respectively. The blended line components 
are indicated by the dotted lines, if any. 
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Fig. 4. — Fitting for the H2 emission lines observed toward HB 21 S. The rest is the same as 
Figure El 
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Fig. 5. — The extinction-corrected population diagrams for HB 21N (left) and HB 21S 
(right). The e xtinctions were corrected, using iV H =iV(H 1)+2N( H.?)= 3.5 x 10 21 cm" 2 
( Lee et all2001) a nd the extinction curve of "Milky Way, Ry = 3.1" ( Weingartner fc Draine 



2001 



Drainell2003l ). The rotational quantum number (J) is printed out near the correspond- 



ing point. 



-22 - 



IC443B 



IC443B 





5.0x10 3 1.0x10" 1.5x10 4 2.0x10" 
Ej (K) 

HB21 N 



5.0x10 3 1.0x10" 1.5x10" 2.0x10" 
Ej (K) 

HB21S 




5.0x10' 



5.0X10 3 1.0x10" 1.5x10" 2.0x10" 
Ej (K) 



Fig. 6. — The model fitting results for the H 2 populations observed in IC 443B (top), HB 
21N (bottom-left) , and HB 21 S (bottom-right) . The IC 443B cases are shown for two different 
parameter settings: X# = —1.7 (top-left) and X# = —2.5 (top-right) (cf. Table HJ). The 
filled-red and open-black symbols indicate the model and data points, respectively. (See 
text for the model description.) The rotational quantum number (J) is printed out near the 
corresponding point. 
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Fig. 7. — The contour of \ 2 i n the plane of model parameters for IC 443B: X# = —1.7 
(top) and X# = —2.5 (bottom) (cf. Table H] and Fig. [6]). The 68% and 95% confidence levels 
are outlined. The tick marks along the contours indicate the directions that \ 2 values are 
decreasing. indicates the model parameter values whose x 2 is minimum; i.e., the best 
fit. 
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Fig. 8. — The model fitting results for the AKARI-IKC spectrum observed in IC 443B. 
The fittings are done for two different parameter settings: X# = —1.7 (top) and X# = —2.5 
(bottom) (cf. Table |5]). The histogram, solid-line, and dotted-line indicate the spectrum data, 
the fitting result, and the adopted continuum, respectively. 
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Fig. 9. — The contour of \ 2 i n the plane of model parameters for IC 443B: X# = —1.7 
(top) and X# = —2.5 (bottom) (cf. Table [5] and Fig. [8]). The 68% and 95% confidence levels 
are outlined. The tick marks along the contours indicate the directions that \ 2 values are 
decreasing. indicates the model parameter values whose x 2 is minimum; i.e., the best 
fit. 
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Table 1. Summary of the AKARI IRC Observations 



Region 


Pointing Position 


Observation ID 


AOT a 




(RA, Dec; J2000) 






HB 21N 


(20:47:26.891, +51:15:04.80) 


1420797. [1-4], 1420798.1 


IRCZ4 


HB 21N BG 


(20:47:32.237,+51:16:34.86) 


1420799. [1-2] 


IRCZ4 


HB 21S 


(20:46:10.398,+49:58:15.20) 


1420800. [1-4], 1420801.1 


IRCZ4 


HB 21S BG 


(20:46:04.220,+49:58:13.61) 


1420802. [1-2] 


IRCZ4 



a Astronomical Observation Template. It is a pre-defined observation sequen ce. 
See the IRC Data User Manual for Post-Helium Mission. (jOnaka et al.ll2009i ). 



Not e. - The observation summary for the IC 443B data is in IShinn et al 

(boiih . 
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Table 2. Observed H 2 Emission Lines 



Transition 


Wavelength 


Observed Intensity 




(H 


(10" 6 erg 


s 1 cm 2 sr r ) 






HB 21N 


HB 21S 


v = 1-0 0(3) 


2.80 


9.1±2.0 


11.6±2.0 


v = 2-1 0(3) 


2.97 


<4.1 




v = 1-0 0(4) 


3.00 


<5.7 




v = 1-0 0(5) 


3.23 


<6.2 


8.3±1.7 


v = l-o 0(7) 


3.81 


<3.8 




v = 0-0 S(13) 


3.85 


<4.9 




v = 0-0 S(ll) 


4.18 


8.4±1.4 


9.0±1.3 


v = 0-0 S(10) 


4.41 


7.5±1.4 


8.7±1.3 


v = 0-0 S(9) 


4.69 


19.7±2.4 


21.6±2.5 



Note. — For those lines whose significance is lower than 
3.0, the intensities are expressed with 90% upper confi- 
dence limits. The error includes both statistical and sys- 
tematic components. See text for detail. 
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Table 3. Extinction Corrected H 2 Column Density 



Level 


Level Energy 


log N(H 


2; v,J) 


(v,J) 


(K) 


(cm 


- 2 ) 






HB 21N 


HB 21S 


(0,11) 


10261. 


15.10±0.05 


15.14±0.05 


(0,12) 


11940. 


14.50±0.08 


14.56±0.07 


(0,13) 


13703. 


14.39±0.07 


14.42±0.06 


(0,15) 


17443. 


<13.90 




(1,1) 


6149. 


14.63±0.10 


14.74±0.07 


(1,2) 


6471. 


<14.62 




(1,3) 


6951. 


<14.82 


14.95±0.09 


(1,5) 


8365. 


<14.97 




(2, 1) 


11789. 


<14.13 





Note. - The extinctions were correc ted, using 
N„=N(R I)+2JV(H 2 )= 3.5 x 10 21 cm" 2 Lee et al. 



200ll ) and the extinction curve of "Milky Way, R 



v 



3.1" f lWeingartner & Draindl200ll ; iDraind 120031 ). 
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Table 4. Level Population Fitting Results for the Model Parameters 



Region 


log[iV(H 2 )] 
(cm -2 ) 


log[n(H 2 )] 
(cm -3 ) 


b 


X 

(= log 


a 

H 

'n(HI)" 
n(H 2 ) 


OPR 

) 


(= X 2 /d.o.f) 


IC443B 
IC443B 
HB21N 
HB21S 


18.6±C;1 

is.8±8:i 

17.4 
17.7 


e o+0.4 
O.8_ 02 

6 5 +0 ' 3 

O.O_ 02 

5.5 
5.5 


o.y_ .4 

4 ^+0.5 
4 -°-0.5 

2.6 
3.3 




1.7 

2.5 
2.0 
2.0 


o i +0.2 
^•-■■-0.2 
9 9+O.2 
z - z -0.3 

1.8 

1.6 


1.8 (=7.2/4.0) 
1.4 (=5.5/4.0) 



a The parameter Xh is fixed as shown. See section 13.21 for more detail. 

Note. - The confidence limits are given with a 68% significance (cf. Fig. [7]). N(H.2) 
means iV(H 2 ; T > 1000 K). See section I3T21 for the detailed description about the parame- 
ters. 
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Table 5. Spectral Fitting Results for the Model Parameters 



Region 


log[iV(H 2 )] 
(cm -2 ) 


log[n(H 2 )] 
(cm -3 ) 


b 


X 

(= log 


a 

H 

' n(HI) " 
n(H 2 ) 


OPR 

) 


A,i/ 

(= X 2 /d.oi) 


IC443B 
IC443B 


1 8 Q+0.02 

i s q +0 - 02 

io - a -0.02 


c 5+o.n 

J - J -0.04 

£ O+0.03 

u -°-0.05 


4 fi +0 - 21 

4 -°-0.05 
c 4+O.O8 
d -^-0.05 




1.7 
2.5 


2 4+0.3 

z -^-0.2 

2 4+0.3 

Z '^-0.2 


1.93 (=408.7/212.0) 
1.89 (=401.4/212.0) 



a The parameter Xh is fixed as shown. See section 13.21 for more detail. 

Note. - - The confidence limits are given with a 68% significance (cf. Fig. [9]). iV(H2) means 
iV(H 2 ; T > 1000 K). See section I3~2l for the detailed description about the parameters. 
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Table 6. H 2 Emission Lines from the Best Fit Models of the IC 443B data 



Transition 


Wavelength 


Modeled Intensity 




(/mi) 


/ -i r\ f\ 1 9 

(10 erg s cm 


sr ) 






X H = -1.7 X H = 


-2.5 


v =2-1 Q(9) 


2.72 


1 7 

X . 1 


1 3 

X . O 


v =1-0 Q(13) 


2.73 


8 


7 


v =2-1 0(2) 


2.79 


2 1 


1 7 

X . 1 


v =1-0 0(3) 


2.80 


uu.y 


rq n 


v =2-1 0(3) 


2.97 


6 Q 

VJ . <J 


5 8 

*J . (J 


v =1-0 0(4) 


3.00 


98 5 


2Q 2 


v =2-1 0(4) 


3.19 


3 1 


2 6 


v =1-0 0(5) 


3.23 


^ Q 


OG.O 


v =2-1 0(5) 


3.44 


6 2 

VJ . 


5 1 

tj . X 


v =1-0 0(6) 


3.50 


15.5 


15.3 


v =2-1 0(6) 


3.72 


1.9 


1.6 


v =1-0 0(7) 


3.81 


21.6 


20.5 


v =0-0 S(13) 


3.85 


15.2 


16.0 


v =0-0 S(12) 


4.00 


11.2 


11.4 


w =2-1 0(7) 


4.05 


2.9 


2.3 


v =1-0 0(8) 


4.16 


4.8 


4.5 


v =0-0 S(ll) 


4.18 


44.6 


45.0 


v =0-0 S(10) 


4.41 


30.5 


30.9 


v =1-1 S(ll) 


4.42 


4.6 


4.1 


v =2-1 0(8) 


4.44 


0.7 


0.6 


v =1-0 0(9) 


4.58 


5.4 


4.8 


v =0-0 S(9) 


4.69 


111.0 


117.6 



Note. — The best fit models are from the spectral fitting 
(Figure E]) and the extinction effect is applied. 
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Table 7. Comparison between the Expected and Observed Intensities of Bra 



Region 


Expected 3, 


Observed b 




(erg s _1 cm" 


2 sr- 1 ) 


IC443B 


(9.3e-07, 5.6e-05) 


<1.7e-05 


IC443C 


(3.1e-06,1.9e-04) 


<1.8e-05 


IC443G 


(1.8e-06, l.le-04) 


<2.0e-05 


HB21N 


(1.3e-07,7.8e-06) 


<4.8e-06 


HB21S 


(1.4e-07,8.4e-06) 


<4.7e-06 



a The expected intensity range of Bra 
is calculated from the obs erved inten- 
sity of v = — S(ll) (jShinn et al. 



20 111 . Table [2]), emp loying the dissociative 



J-shock model of Hollenbach fc McKee 



(119891 ) . The values are calculated for no 



10 s 



cm 



and v , = 40 — 80 km s 



b 99 % upper confidence limit. 



